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Abstract: Using the method of prolongation we generate new stationary solutions from a 
recently obtained simple particular solution for relativistic transverse flow with cylindrical sym- 
metry in 1+3 dimension. This is an extension of the longitudinal Bjorken flow ansatz and can 
be applied to situations during a first order phase transition in a fast expanding system. The 
prolongated solution allows us to trace back the flow profile from any transverse flow conjectured 
^1 ■ at the end of the phase transition. 
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Introduction 

Hydrodynamics often allow for nonrelativistic scaling solutions. Relativistic flow, how- 
ever, seemed long to be an exception: besides Bjorken's 1+1 dimensional ansatz and the 
spherically symmetric relativistic expansion, no analytical solution was known [§]. 



In a recent paper we presented an extension of Bjorken's ansatz || for longitudinally 
and transversally relativistic flow patterns with cylindrical symmetry in 1+3 dimensions. 
It satisfies u^d^u u = as well as the original Bjorken ansatz. This analytical solution of 
the flow equations of a perfect fluid is in particular valid for physical situations when the 
sound velocity is zero, c 2 s = dp/de = 0, with energy density e and pressure p(e). 

In particular this happens during a first order phase transition, the pressure is constant 
while the energy density changes (in heavy ion collisons increases and drops again). This 
should, in principle, be signalled by a vanishing sound velocity. A remnant of this effect 
in finite size, finite time transitions might be a softest point of the equation of state, 
where c 2 s is minimal. In fact, this has been suggested as a signal of phase transition by 
Shuryak[|J, and investigated numerically in several recent works || ||. 

In this paper we generalize that simple analytic solution further by exploring the symme- 
tries of the nonlinear partial differential equation determining the flow at the softest point 
(the relativistic Euler's equation). The method of prolongation ensures us that all possi- 
ble transformations of dependent and independent variables can be found which comply 
with the equations to be solved by solving linear partial differential equations only. Then 
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from any simple (or trivial) particular solution of the original nonlinear problem several 
new classes of solutions can be generated. In this article we present the simple solution 
published recently, we review the method of prolongation and obtain a general class of 
transverse flow patterns for the relativistic problem at the softest point. We use this 
solution for calculating backwards from a conjectured final state at CERN SPS ||. 

Nonrelativistic analytic solution has been also given several times, with respect to heavy 
ions see [|,[K|[n|. 



Simple relativistic transverse flow at the softest point 

The 1+1 dimensional Bjorken flow four-velocity is a normalized, timelike vector. It is 
natural to choose this as the first of our comoving frame basis vectors (vierbein). The 
three further, spacelike vectors will be constructed orthogonal to this, separating the 
two transverse directions. This basis fits excellently to a cylindrical symmetry and to 
longitudinally extreme relativistic flow. 

<£ = (7.7.0.°), 4 = (7. 7.°.°), 

^ = (0.0.7.7.). »J = (0.0.-7.7.). (1) 

with t time coordinate and z longitudinal (beam-along) coordinate, x and y transverse, 
cartesian coordinates. The cylindrical radius is given by, r = yjx 2 + y 2 , and r is the 
longitudinal proper time: r = y/t 2 — z 2 . 

We consider ideal fluids ("dry water"), where the energy momentum tensor is given by 

T^u = (e + p)UfjU u - pg^, (2) 

and the equation of state is given in the form of p(e). The ansatz for an almost boost 
invariant flow with some transverse, cylindrically symmetric component is then given by 

u» = l(el + vel), (3) 

using the Lorentz factor 7 = (1 — v 2 )^ 1 ^ 2 . This four- velocity is normalized to one: 

= 7 2 - 7V = 1. (4) 

In a recent publication we formulated Euler's equation in terms of this ansatz and coordi- 
nates. For a situation with zero comoving gradient of the pressure it has been simplified 
to a sole partial differential equation for the transverse flow velocity component v (r, r): 



This equation is valid even for relativistic transverse flow, it is a quasilinear partial dif- 
ferential equation. The separable solution,?; (r, r) = a(r)b(r), which is regular at the 
cylindrical axis r = is given by 

v= r -. (6) 

In a recent paper we then explored consequences of such a linear transverse velocity profile 
for the time evolution of the energy density at constant pressure, i.e. the quark matter 
to hadron matter phase conversion pi]. 



Our concern now is to construct a more general solution of eq. ([5|) , which is able to fit any 
flow profile at a given longitudinal time 

v(r ,r) = v {r). (7) 

In order to explore the quark matter hadron matter phase conversion we utilize simple 
bag model equation of state, as was done in [|lj. The energy density of x P ar t quark 
matter and (1 — x) hadron matter is then given by 

e=(a q T* + B) X + a h T*(l-x), (8) 

at the phase equilibrium temperature T c , determined from the equity of quark and hadron 
pressure. The resulting evolution equation turns to be 

+ udX = 0, (9) 

with 

» = do) 

Here the ratio of Stefan-Boltzmann constants a Qi h depends only on the relative number 
of the effective degrees of freedom. By comparing a relativistic pion gas to quark gluon 
plasma ah = 3 and o q = 37 and one obtains v w 0.12. 



Now using the flow ansatz discussed above and eq.(|^) one arrives at 

( d d\ , . f dv v 1 \ 

(* + "*)* + <* + % + ; + 7]=°- < u > 

Furthermore, expressing the quark matter part as 

X{r,r) = -U+—Z{r,r), (12) 

rr 

we obtain a simpler equation for Z, 

£ z+ lw = ' < 13 > 

which formally resembles a one-dimensional conservation equation with density Z and 
convective current vZ. It is noteworthy that for any solution Z of this equation (|T3| ) 



3 



f(v)Z is also a solution with f(v) being an arbitrary (but differentiable) function of the 
transverse velocity v. This freedom may be used for fitting a quark matter percentage 
profile at the same time when the transverse velocity profile is known, especially the x — 
(end of phase conversion) situation may be traced back up to a time when somewhere 
X = 1 (full quark matter) is reached. 



Eq . (|13D is formally a conservation equation, whose solution, Z(r,r) satisfies 

Q(r) = [ R(T) Z(r,r)dr = Q(r ). (14) 
Jo 

The longitudinal time derivative of this expression vanishes leading to 

dO ■ r R ( T ) dZ 

*=J*r)Z(r,r) +J ( §^ = 0. (15) 

Replacing dZ/dr from eq. (|l3D we arrive at 

'R(t)-v(t,R(t))]z(t,R(t)) = (16) 

if we take into account the regularity condition v(t, 0) = 0. This is an evolution equation 
for the transverse radius R(t) of the phase mixture region where Z ^ 0, 

R(t) = v(t,R(t)). (17) 

It is a first order ordinary differential equation which can be easily solved once v(r,r) is 
known. 

The conserved quantity Q(r) can be related to the spatial average of the phase mixture 



ratio by using eq . ([12|) : 

2 r R 

(x) = -jp f_ x(r,r)rdr. 



o 



Finally we obtain 

{X) = - U+ Vlkr) m - (19) 



In the followings we turn to a brief presentation of the method of prolongation, which 
shall be used for generalizing the particular relativistic transverse flow solution presented 
in ill. 



The method of prolongation 



The method of prolongation "prolonges" symmetries of a (system of) nonlinear partial 
differential equations, to transformations of dependent and independent variables and all 
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partial derivatives of the dependent variables up to a degree less then the rank of the 
equation itself. For an incomplete list of references about the application of this method 
in physics see Jl2[ . 



Let us generally denote the set of independent variables by x, the dependent ones by v, 
partial derivatives by v x , v xx etc. The equation(s) to be solved let be given by 

A{x,v,v x ,...) =0. (20) 

Transformations on x and v (generating corresponding transformations on the partial 
derivatives) may let this equation remain valid in the transformed variables as well. These 
transformations are symmetries of the equation, and they can be viewed as integrals of 
infinitesimal transformations generated by the vector field (Killing vectors) 

K = $(u,x)-^ + Z(v,x)-^. (21) 



This relation can be easily understood by inspecting infinitesimal transformations, 

v = V + E&, 



x = x + e£. (22) 



In this case the derivative transforms to 



dv i _ d<E> 

^+^[^-^)+0(e 2 ) (23) 



dv _f x +ef x dv , Jd<$> d^dv 



dx 1 + e^f- dx \ dx dx dx 

dx v 

The total derivatives are 

ax 

£ = ZvVx + Zx. (24) 
The differential equation, A = can be expanded around the original solution 

A(v,x,v x ) =A(v,x,v x )+e$—A + et—A + el — --±v x \ —A + <D(e 2 ). (25) 

The second and third term of O(e) in the above expansion resembles KA, but all terms 
to this order constitute a Killing field in the extended (prolongated) space spanned by 
v,x and v x . This is the (first) prolongation of the symmetry 

d 

W {1) K = K + § x — (26) 

dv x 
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with 

d$ d£ d _ , i . „. 

$ ^ - u *d! = & ($ - v *° + (27) 

Since both u(:r) and -u(x) are solutions, A(t>, x, v x ) = and A(v,x,v x ) = for arbitrary 
e and we arrive at 

pr«K(A)| A=Q = 0. (28) 



The requirement of transforming only among solutions leads to a number of linear partial 
differential equations for the unknown functions $ and £ (which may contain, in general, 
several components) by equating the coefficients of each monomials in v, v x etc. with zero. 
Eq . (|28|) is a property of the transformation and therefore is fulfilled for any solution of 
the original equation. The A = constraint is important, it reduces the dimensionality of 
the space of the coefficient functions, usually by expressing one of the highest oder partial 
derivatives in terms of others. The resulting system of equations in principle can be solved 
(being linear). The corresponding vector field K will then contain in general a number 
of constants or undetermined functions selecting out classes of transformations which 
generate new solutions from a given particular solution (BO). The method of prolongation 
ensures us that we have considered all possible transformations. 



New solutions 

Now we rewrite equation (|5|) with notations of the method of prolongation: 

A = v T + vv r = 0, 

Odd 
K = $(v,T,r)— + T(v,T,r)— + R(v,T,r)—, 

pr {1) K = K + $ T — + <F— . (29) 

Applying the prolonged Killing vector field to the nonlinear equation to be solved we get 

pr (1) KA = $ r + v<$> r + v r $ = 0. (30) 
The prolonged coefficients of the Killing vector field are 

$ r = ($ T + v T $ v ) - (T T v T + T v v 2 T ) - (R T v r + R v v T u r ), 

$ r = ($ r + v r $ v ) - (T r v T + T v v T u r ) - (R r v r + R v v 2 r ). (31) 
Replacing it to the prolongation condition eq. (|30|) , we arrive at 

pi (1) K (A)| A=0 = $ T + v$ r + tv($ - R T ) + vv r (T T - R r ) + v 2 v r T r = 0. (32) 

6 



This has to be fulfilled for arbitrary v(r,r). In this equation, since we have already used 
the A = constraint (the original equation to be solved) for eliminating v T = —vv r , each 
coefficient in the polynomial of v and v r vanishes separately, leading to a number of linear 
partial differential equations: 

$ r = 0, $ r = 0, $ = R T , T T = R r , T r = 0. (33) 

To these equations (being linear) a general solution can be given: 

T(v,T,r) = b(v) + a(v)r, 

R(v,r,r) = c(v ) + a(v )r + <f>(v)r, 

<5>(v,r,r) = 4>(v). (34) 

There are four undetermined functions in this general solution to the Killing vector prob- 
lem, hence there are four general classes of transformations leaving the fulfillment of the 
original equation untouched. The corresponding vector fileds, which generate infinitesimal 
transformations, are 

K ' - «">(! + '!)- 

A'3 = 

Kt = c(v)^-. (35) 

A general finite transformation of any solution of A = is given by 

(v, f, f) = exp {eiKx + e 2 K 2 + e 3 K 3 + e 4 K A ) (v, r, r). (36) 

Three of these transformations are quite trivial: K± generates translation of the variable 
r by c(v ), K 3 generates translation of the variable r by b(v), and K 2 a strecth of both r 
and r by e a( - v \ Only K\ involves a transformation of v as well: 

v = v + ip(v), 

f = T, 

f = r + rtp(v), (37) 

with 

dv 



™ = fwr (38) 
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Combination of these transformations generates from the simple particular solution v = 
r/r the following (implicit) solution 

v - r + A{V) (39) 
Fitting this form to an initial transverse velocity profile, leads to 

„ 0= „ (0 , r)= r +' 4( ;° ( : )) - <«> 

B(v (r)) 

which is equivalent to 

r = v B(v )-A(v ) = f(v ). (41) 
Finally we arrive at the general solution 

r-vr = f(v), (42) 

with f(v) being the inverse function of the initial transverse velocity profile vo(r). The 
general solution is implicitely given by 

v — Vo(r — vt). (43) 

It is easy to verify that this implicit form still satisfies the original equation (13). 



Finally let us consider some possible initial transverse velocity profiles. For the linear 

case, 

t>o(r) = ar, (44) 

we obtain 

v = a(r — vt) (45) 
which can be resolved to the explicit expression 

ar r 



1 + ar r+l/a 



(46) 



presented already in |lj . Initial profiles with nonlinear or non-invertable functions can be 
resolved only numerically, but a solution of implicit equations are usually faster and more 
stable than algorithms for solving partial differential equations numerically. 

Flow at SPS 



A particularly interesting enterprise is to start with a transverse flow profile conjectured 
from experimental data and evolve the flow pattern backwards in (longitudinal) time. For 
the CERN SPS heavy ion experiments Ster, Csorgo and Lorstad have proposed recently 
such a profile based on one-particle transverse momentum spectra and HBT correlation 
measurments [f|. The reconstructed flow in the notations of the present paper can be 
written as, 

jv = sinhr/t = (u t )-f- = ar, (47) 
Kg 
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Figure 1: Relativistic transverse velocity profiles from the prolonged analytic solution. 
The curves are results of a time backward calculation from a conjectured final state at 
t = 6 fm/c with a radius of R = 7.1 fm. The average portion of the quark matter in the 
mixed phase is indicated at the respective values of 0, 1/3 and 1. 

at a late time instant r = 6 fm/c. The fit by Ster et. al. concludes at (u t ) = 0.55 ± 0.06 
and Rq = 7.1 ± 0.2 fm. This gives rise to the following proportionality coefficient: a = 
0.0775 ± 0.01. This suggestion has improved on the linear transverse flow assumption, it 
nowhere exceeds the speed of light. 

Fig.l shows the reconstructed relativistic transverse velocity profiles ending at the one 
suggested by Ster et.al. It was assumed that in the final stage the phase transition just 
ended in an overall hadronic phase. 

As insertions show the mixed phase should have started then as early as r = 1.14 fm/c, 
but the quark matter part was decreased to 1/3 on the average soon (r = 2.4 fm/c). 
The velocity profiles do not differ from linear qualitatively in the range of interest. It is 
more pronounced how the radius of the mixed phase expands from R$ = 5.4 fm to the 
reconstructed final stage with R = 1.1 fm. 



Conclusion 

With the help of the method of prolongation we explored the symmetries of Euler's 
equation for a relativistic transverse flow at the softest point (zero comoving gradient 
for pressure). We arrived at an implicit solution (|13| ) which presents a correspondence 
between the initial transverse velocity profile and that at arbitrary longitudinal proper 
time. Although in a general case an implicit equation can only be solved numerically, it 
is a more elegant, stable and concise way to solve the flow problem than the numerical 
integration of partial differential equations. 
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With this perspective a transverse flow pattern obtained from measurements in heavy ion 
experiments at CERN SPS, where probably the mixed quark-gluon and hadronic phase 
has been realized, can also be re-calculated at the beginning of phase conversion (pure 
quark matter). In turn this might help for obtaining an improved insight into further 
characteristics (temperature, effective masses) of the quark matter at CERN. 
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